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ince the invention of rotating 
a machines, the pursuit of higher 
power output has driven machine 
speeds higher and higher. With the 
hreaking of the first balance resonance 
“barrier” {achieved by DeLaval with a 
steam turbine in 1895), rotating 
machines were shown Lo be able to oper- 
ale continuously above the first balance 
resonance. However, with this new capa- 
bility came a new problem for machines 
using fuid-lubricated journal bearings: 
Muid-induced instability. Now recog- 
nized nol only in bearings, but in seals 
and in the main flow of fHluid-handling 
machines, fluid-induced instability has 
been a major focus of research efforts 
around the world as engineers have 
struggled to achieve higher operating 
speeds and higher power, while at the 
sume time maintaining rotor stability, 
Over Lhe years, many different methods 
and approaches have been developed by 
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researchers to identify and under- 
stand the important parameters 
that influence rotor stability. 

In Walter Evans’ book, 
“Controal-System Dynamics,” |Ref. 
3| is a powerful tool for the visual- 
Imation and analysis of rotor stabil- 
ily: the Root Locus. Don Bently 
worked with Walter Evans at the 
time that this hook was first pub- 
lished. The Root Locus technique 
i$ a powerful, graphically-based method 
for presenting rotor stability informa- 
tion, and, with the arrival of more power- 


ful personal computers, is becoming 


Increasingly casy to implement. The 
Rool Locus is based on very simple con- 
cepts and can be easily understood and 
used for visualizing rotor system insta- 





bility and parameter changes to 
enhance the stability margin. 


The Root Locus 

Use of the Root Locus technique 
requires a mathematical model of the 
rotor system under study. Although 
complicated, multi-Degrce-of- 
Freedom models can be useful for the 
study of particular rotor system stability 
problems, a simple rotor model can pro- 
vide a great deal of insight into the 
behavior of more complicated systems. 
The simple model presented in [Ref.5] 
Will be used here to initially demon- 
strate the Root Locus procedure. The 
model leads ta the 2nd-order, linear, 
differential equation, 


Mz+D2+(KjDasiz=0 (D 





Figure 1 


s-plane graph of the root loci of the two roots of the simple rotor model, These plots were 
made by varying rotor rotative speed, 9, from 0 to 2000 rac)'sec. The zero speed points are 
marked with an “X.” The vertical axis marks the Threshold of Stability, The left half-plane 
corresponds to stable roots; the right half-plane corresponds to unstable roots for the 
system under consideration, The lower root poses no threat to stability, but the upper root 
crosses the vertical axis and moves to the right half-plane of instability. 
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where z = xt+yy is the complex 
displacement 


coordinate, 

Ml = rotor system modal 
mäss, 

D = rotor system modal 
damping, 

K = rotor system modal 
stiffness, 

ÀA = [lwd average circum- 
[erential velocity ratio, 

2 = rotor rotative speed, 

j = -1 


and the dots indicate derivatives with 
respect tọ time, £. Mote thatx and y are 
the displacement coordinates in the 
Plane perpendicular to the rotor axis of 
rotation, 

In Rool Locus analysis, the external 
forcing of the system is assumed to be 
Zero in order to allow examination of 
the free vibration characteristics of the 
system, A standard solution technique 
for differential equations yields the 
Characteristic Egteation of le system, 


Mis* + Ds + K-jDAn=0 (2) 


Because Equation (2) 1s quadratic, 
there will be two roots. Each of the two 
roots of this equation are complex 
numbers af the form, 


SS + jaw (3) 


where œ and w are complicated fung- 
tions of, DA , K, and &. 


Verlical axis crossing mak 
the Threshold of Instability 
LHH & 1050 radrsec 
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Figure 2 


Important stability root of Figure 1. Rotative speeds are shown in 
rad‘sec. The Threshold of Stability is shown to be near 1050 rad, 


SEC. 





Each of the two roots has a direct and 
quadrature part. (While the direct and 
quadrature parts of the roots corre- 
spond to the “Real” and “Imaginary” 
parts of the roots as typically described 
in mathematical theory, there is nothing 
“Imaginary about the meaning of these 
equations.) The two roots lead to a geti- 
eral solution. bor purposes of illustra- 
tion, using only one root, the behavior of 
the rotor model can be expressed as 


z= Zeateivl (4) 


where Z is a complex rotor position con- 
slant which is nol important in the Root 
Locus stability analysis. 


The rotor natural frequency (the 
angular velocity of precession) is 
defined by w, which, while a function of 
the rotative speed, @, Is usually not 
equal to the rotative speed. 

The rate of growth or decay of the 
amplitude of the orbit of precession, z, 
is controlled by œ. It controls the decay 
rate or growth rate of the rotor orbit 
indirectly, through the first exponential 
function in Equation (4) and thus con- 
trols the transient response amplitude 
of the rotor system due to some distur- 
bance. If a<0, the first exponential 
function will become smaller with Lime, 
and the precession orbit will decay with 
time back to the original equilibrium 
postion. However, if a>0, the preces- 
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sion orbit will increase in size in an 
unbounded fashion! 

The roots of the characteristic egua- 
lion completely define the two time- 
dependent characteristics, growth (or 
decay) and precession frequency, of the 
free vibration of the rotor system. The 
Iwo Components, g and w, of cach rool, 
s, can be graphed in the s-plane. By 
systematically varying a single rotor 
mode] parameter, curves in the s-plane 
can be generated, Rach curve is a locus 
of roots of the characteristic equation, 
or a Root Locus curve and 1s a function 
of the chosen parameter. 


Unbounded growth of the rotor orbit 
is the condition of rotor instability, 
Thus, for stability, all a<0; that is, all 
roms of the characteristic equation musi 
he ta the felt half of the s-plane. The con- 
dition when any a= defines the 
Threshold of Stability. 


Such a pair of curves for the two roots 
of the rotor model are shown in Figure 
l. These Root Locus curves were gener- 
ated by varying the rotative speed, Q, 
from Ü to 2000 rad/sec (in the counter- 
clockwise or positive angular direction}, 
The large “X” on the plots marks the 
zero rotative speed point. The roots are 
plotted with the Growth or Decay 
parameter, a, as the horizontal axis, and 
the natural frequency or precession 
rate, w, on the vertical axis. Note that 
one of the roots is above the horizontal 
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Example of a Root Locus analysis of a rotor system different than 
the one presented in Figures 1 and 2. Increasing the direct stiff- 


ness, K, from 4000 Ibin to 8000 lbn results in an increase in the 
Threshold of Stability of more than 25%, 
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axis, indicating forward precession, 
while the other rool is below the hori- 
zontal axis, indicating reverse preces- 
sian, Note also that the reverse 
precession root moves toward greater 
decay rate with increasing rolative 
speed while the forward precession rool 
moves toward lower decay rate (lower 
absolute value of a) with increasing 
rolalive speed, . 

Most importantly, the forward root 
does cross the vertical axis at some 
speed. This indicates that, as a becomes 
equal to and then greater than zero, the 
rotor system will became unstable and 
proceed into a new limit cycle of either 
whirl or whip. In this model, the reverse 
root will never pose a threat to stability 
for positive, counterclockwise rotor 
rotation, it would pose a threat to stabil- 
ity for negative, clockwise rolor rota- 
tion, For that reason, ef! roots of a rotor 
model should be examined for potential 
ZETO CrOSSINES. 


A graph showing only the important 
stability root of Figure | is presented in 
Figure 2. In this figure, the numbers 
along the curve represent rotor rotative 
apeeds in rad/sec, The Threshold of Sta- 
bility is shown to be approximately 1050 
Taclsec, 


Anexample ofa Root Locus analysis 
is shown in Figure 3, In this example, a 
different rotor model than the one pre- 
sented in Figures 1 and 2 was used. The 
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rotor rotative speed was varied from 0 to 
400 rad/sec, The Threshold of Stability 
lor this rotor system was found to be 
about 175 rad/’sec. At that point, the 
rotative speed was fixed at 180 rad/sec 
and the direct stiffness, K, was varied 
from its original value of 4000 Ib/in to a 
inal value of 10000 Ibin. The improve- 
ment in rotor stability can be clearly 
seen as the Root Locus crosses back into 
the left-hand side (the stable side) of 
the 4-plane. 

Note also that the rotor natural fre- 
quency (precession frequency) 
Increases, as would be expected from an 
Increase in stiffness of the system. 
Finally, the direct stiffness was fixed at 
BLU Ibin, and the rotor rotative speed 
was increased from 180 rad/sec to 300 
rad/sec. The new Threshold of Stability 
was found to be about 220 rad/sec, a 
25% increase. This could reflect the 
improvement due to increasing bearing 
fluid pressure or io a deliberate, 
“friendly” misalignment placing the 
shall into a higher eccentricity operat- 
Ing position. 

A tinal example is a Root Locus anal- 
ysis Of a two-stage compressor. The 
compressor rotor was supported by twa 
Nuid-lubricated bearings and had a seal 
located at the mid-span between the 
two stages. The compressor had a his- 
tory of whirling under certain startup 





conditions at a speed between the first 
ind second lateral balance resonances. 
It was decided to compare the relative 
effectiveness of anti-swirl injection (see 
Ref, |10]} at the center seal to anti-swirl 
injection at the fluid-lubricated 
bearing. 

A more complicated model was 
developed for this problem which 
resulted in a 6th-order Characteristic 
Equation with complex coefficients, 
The roots of this equation were found 
numerically using a Siljak polynomial 
method, 

Since the characteristic equation was 
Oth-order, there were six roots. An 
s-plane Root Locus plot of the four 
rightmost roots of the model ts pre- 
sented in Figure 4. The other two roots 
were located far to the left in the s- 
plane and had no importance in the 
stability analysis. To create the Root 
Locus plots for this model, the shaft 
rotative speed was varied from (to 4000 
rad/sec. Note that one of the roots 
crosses the vertical (precession) axis 
into the right, unstable half-plane. The 
crossing point defines the Threshold of 
Stability for this model. An enlarged 
view Of this important stability root is 
shown in Figure 5. The rotor system is 
predicted to become unstable at a rota- 
live speed of about 1050 rad/sec. 

To examine the effectiveness of anti- 
swirl at the bearing, Ap, the center seal 
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Figure 3 


Root Locus s-plane plot of four of the six roots of the compressor 
model for rotative speed fram 0 (marked with an “X™ to 4000 rad; 
sec. The other two roots are located far to the left and do not have 
any importance in the stability analysis, One of the roots crosses 
the vertical axis into the right, unstable half-plane. The other roots 
pose no threat to stability in this speed range. 
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Enlarged view of the important stability root of Figure 4. The 
rotative speeds in rad/sec are indicated on the curve. The Thresh- 
old of Stability is approximately 1050 radsec. Note the similarity 
of this Root Locus to the one in Figure 2 which was generated 
using similar rotor parameters using the model of Equation (1). 
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average circumferential velocity ratio, 
was held constant at 0.48 while Ag, the 
Tluid-lubricated bearing average cir- 
cuntlerential velocity ratio, was varied 
from 0.5 100.1 in five steps. To examine 
the effectiveness of anti-swirl at the cen- 
ter seal, Ag was held constant al 0.4% 
while As was varied from 0.5 to 0.1 in 
live sleps, 

Figures 6 and 7 show the effect of the 
two anti-swirl injection strategies on the 
Important stability root of Figure 3. 
These Root Locus plots immediately 
reveal that anti-swirl is much more 
effective when applied to the center seal 
than when applied to the fluid- 
lubricated bearing. Even a small 
amount of anti-swirl injection at the 
seal resulis in a substantial Improve- 
ment in the Threshold of Stability. 


Conclusion 

[tis important to point out that the 
Threshold of Stability predicted by the 
Root Locus (and by any method of anal- 
ysis, for that matter} is only as good as 
the model used in the analysis. Because 
of this, any stability analysis based on a 
model of a complex machine is unlikely 
lo predict the exact Threshold of 
Stability. 

The power of the Root Locus method 
comes from its ability to reveal how sen- 
sitive a rotor instability problem is to 
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Figure 6 


Effect of anti-swirl injection at the fuid-lubricated bearing on the 
stability of the compressor. For five values of Ap (0.5 to Ô. 1), Ac 
was held constant at 0.48, and rotor rotalive speed was varied 


from Ô te 2000 rad‘sec. 
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different rotor parameters and from its 
powerful visual presentation of those 
relationships. It allows one to consider 
various solution alternatives and look 
lor Lhe most cost-effective solution to a 
potential rotor stability problem. 
Although complicated, multi- 
Depree-ol-Preedom models can be use- 
ful for the study of particular rotor sys- 
tem stability problems (such as the 
compressor problem above), a simple 


rotor model can provide a great deal of 


insight into the behavier of more com- 
plicated systems, Comparison of Fieure 
5 (based on a complicated model) to 
Figure 2 (based on a simple model using 
similar rotor parameters) shows 
remarkable agreement in the predicted 
Threshold of Stability. 

With the powerful symbolic and 
numeric desktop computation tools 
avilable today, the application of the 
Root Locus technique to rotor stability 
analysis has become much eusier to 
implement. The power of this method 
makes it an excellent tool for engineers 
to use at both the design stage and for the 
exploration of potential fixes for stability 
problems encountered in the field. I 


References 
L Bently, D. E., and Hatch, ©. T., “Application of 


the Root Locus Technique to the Analysis of 


Rotor Stability Problems.” Bently Rotor 
Dynamics Research Corporation Report 
GSR, Gently Rotor Dynamics Research Cor- 
poration, Minden, NY, 


a5 
Hao f- 
525 | 
20 
51E 


510 f 


Procession Rata [radisac] 


ES F 


rad'sec. 


b- 


“Anti-swirl at 
‘Center Seal 


. Bently, D. E., 


a Chestnut, F, and Mayer, RAW, “Semvurieche- 


many and Remulaiiny Svsiem Design” John 
Wiley & Sons, incs New York, MY, 1931. 


. Evans, WR. “Coniroihinen Dynamics 


MeGraw-Hil! Book Company, Jne., New York, 
NY, 1954. 


-Kug B.C., "Automatic Conto! Sysiems? Fifth 


Ed. Prentice-Hall. Ine. Englewood Cliffs. NJ, 
ode, 987. 


- Muszynska, AL “One Lateral Mode Isotropic 


Rotor Response To Nonsynchranaus Excita- 
tion. Yon Karman Institute for Fluid 
Dynamics, course on “Rotor Dynamics ard 
kimananin Turbhonachingy, 1992, 


- Muszynaka, A. “Whirl And Whip Rotor’ 


Bearing Stability Problems? Journal OF Sound 
and Vibration, Y. 110, N. 3, 1986, pp 444462. 
and Musaynska, An “Role of 
Circumterential Flow in the Stability of Fluid- 
Handling Machine Rotors” Texas A&M Pilth 
Workshop on Rotonlynamics Instability Prob- 
lems in High Performance Turbomachiners, 
WASA CP 3026, 1988, pp 415-440. 


» Musaviska, A. “Mullimode Whirl and Whip 


in Rotor{Bearing Systems.” Dynamics of 
Kotating Machinery: Proceedings of the See- 
and International Sympoasium an Transport 
Phenomena, Dynamics, and Design of Rota- 
tng Machinery, Part I1, Hemisphere Publish- 
Ing Company. Mew York, MY. 1990, 


» Bently, D. E. and Muszynska, A., “Perturba- 


Hon Study of a Rotor Besring System: bdenti- 
fication of the Oil Whirl and Oil Whip 
Resonances” ASME Design Engineering 
Division Conference and Exhibit on Wechani- 
cal Vibration and Noise, Cinginnati, Obie, 
98S. Paper 65-DE 1-142, 


| Bently, D. E., and Muszynska, ^., “Amt-Swirl 


Arrangements Prevent Rotor Seal Instability’ 
ASME Journal of Vibration, Acoustics, Stress, 
and Reliability in Design, Y. 111. N. 2, 1989, 
pp žole. 





neal 0 5} 100 
Dacar Rate {1S0 


Figure 7 


Effect of anti-swirl injection at the center seal on the stability of the 
compressor, For five values of Ac (0.5 to 0.1), A, was held con- 
stant at 0.46, and rotor rotative speed was varied from 0 to 2000 
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